Observational constraints on Galileon cosmology 
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We study the cosmology of a covariant Galileon field (j> with five covariant Lagrangians and 
confront this theory with the most recent cosmological probes: the type la supernovae data (Con- 
stitution and Union2 sets), cosmic microwave background (WMAP7) and the baryon acoustic oscil- 
lations (SDSS7). In the Galileon cosmology with a late-time de Sitter attractor, there is a tracker 
that attracts solutions with different initial conditions to a common trajectory. Including the cosmic 
curvature K, we place observational constraints on two distinct cases: (i) the tracker, and (ii) the 
generic solutions to the equations of motion. We find that the tracker solution can be consistent 
with the individual observational data, but it is disfavored by the combined data analysis. The 
generic solutions fare quite well when a non-zero curvature parameter tljp is taken into account, 
but the Akaike and Bayesian information criteria show that they are not particularly favored over 
the ACDM model. 
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I. INTRODUCTION 

The quest to distinguish between a cosmological con- 
stant and dynamical dark energy models has been one of 
the main topics in cosmology. The equation of state wde 
of dark energy can be constrained not only by the Su- 
pernovae la (SN la) data [l[ but by the observations of 
Cosmic Microwave Background (CMB) [2[ and Baryon 
Acoustic Oscillations (BAO) d, For constant w DE 
the combined data analysis of CMB+BAO-i-SN la by 
the WMAP group has given a tight constraint u>de = 
-0.980 ± 0.053 (68 % confidence level) in the flat Uni- 
verse However, the present observations still allow 
a large variation of the dark energy equation of state in 
terms of the redshift z. Moreover, inclusion of the cosmic 
curvature further weakens the constraints on wde(z). 

Over the past decade, many dynamical dark energy 
models have been proposed as an alternative to the cos- 
mological constant (see Refs. @ for review). They are 
broadly classified into two classes-(i) Modified matter 
models, and (ii) Modified gravity models. In the class (i) 
the accelerated expansion of the Universe is induced by a 
modified matter source, whereas in the class (ii) the large- 
distance modification of gravity gives rise to the cosmic 
acceleration. The representative model of the class (i) is 
quintessence based on a minimally coupled scalar field Q , 
but in general there is a degeneracy around u>de = — 1 
if we constrain the quintessence potential from observa- 
tions. 

The modified gravity models proposed so far consist 
of f(R) gravity scalar-tensor theory M, the Dvali- 
Gabadadze-Porrati (DGP) braneworld [T(| model, the 
Gauss-Bonnet gravity [ll|, f(R,Q) gravity [l2[, and so 



on. In general we need to recover the General Relativis- 
tic behavior in the region of high density for the consis- 
tency with solar-system experiments (l3j , while the large- 
distance modification of gravity leads to the cosmic accel- 
eration today. Moreover we require that the models are 
free from ghosts and instabilities 
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For example, the 

DGP model is plagued by the ghost problem [lj| in addi- 
tion to the incompatibility with observational constraints 
[lf|. The dark energy models in which the Lagrangian 
includes a general function / of the Gauss-Bonnet term Q 
also result in violent instabilities for small-scale perturba- 
tions [IT}- In f{R) gravity and scalar-tensor theory, the 
functions f(R) or the field potentials need to be carefully 
designed to satisfy the above-mentioned demands [lH . 

In the DGP model a branc-bcnding mode </> (i.e. lon- 
gitudinal graviton) gives rise to a field self-interaction 
of the form U\4>(d 11 4>d a 4>) through the mixing with the 
transverse graviton This allows the decoupling of 

the field <p from gravitational dynamics in the local re- 
gion by the so-called Vainshtein mechanism [2(|. Then 
the General Relativistic behavior can be recovered within 
a radius larger than the solar-system scale. The self- 
interaction □0(<9 A1 0(9^0) satisfies the Galilean symmetry 
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bfj_ in the flat space-time. While the DGP 
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model suffers from the ghost problem, the extension of 
the field self-interaction to more general forms satisfying 
the Galilean symmetry may allow us to avoid the appear- 
ance of ghosts. 

Nicolis et al. [2l| showed that there are only five field 
Lagrangians (i = 1, • • • ,5) that respect the Galilean 
symmetry in the Minkowski background. These terms 
lead to only the second-order field equations and hence 
we do not need to worry about extra degrees of freedom 
coming from higher-order derivatives. 

If the Lagrangians are varied in the curved space- 
time, the terms £4,5 give the field equations higher than 
the second-order. Deffayet et al. [22[ derived the co- 
variant Lagrangians d (i = 1, • • ■ ,5) that result in only 
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the second-order equations, while recovering the Galilean 
symmetry in the Minkowski space-time. This can be 
achieved by introducing field-derivative couplings with 
the Ricci scalar R and the Einstein tensor G vp in the ex- 
pression of £4,5. Since the existence of those terms affects 
the effective gravitational coupling, the Galileon gravity 
based on the covariant Lagrangians (i = 1, • • ■ , 5) can 
be classified as one of modified gravitational theories. 

The full cosmological dynamics including the terms up 
to £5 have been studied by two of the present authors in 
the flat Friedmann-Lemaitre-Robertson- Walker (FLRW) 
background [H, Hi] (see also Refs. [25|- [H] for related 
works). While the field is nearly frozen during the early 
epoch through the cosmological Vainshtein mechanism, it 
begins to evolve at late times to lead to the acceleration 
of the Universe. The stable dS solution can be realized 
by a constant field velocity. 

Refs. [23], [24[ have shown that, for the covariant 
Galileon theory having dS attractors, cosmological solu- 
tions with different initial conditions converge to a com- 
mon trajectory- a tracker solution. Moreover the back- 
ground cosmological dynamics along the tracker can be 
known analytically in terms of the redshift z. The dark 
energy equation of state exhibits the peculiar phantom- 
like evolution: wue = —7/3 (radiation era), »de = —2 
(matter era), and u>de = —1 (dS era). Note that 
this does not imply the appearance of ghosts. In fact 
the viable model parameter space has been found in 
Refs. [23|, [24| from the conditions to avoid ghosts and 
Laplacian instabilities of scalar and tensor perturbations. 

In this paper we place observational constraints on the 
covariant Galileon gravity using the observational data 
of SN la, the CMB shift parameters, and BAO. In par- 
ticular we derive a convenient analytic formula for the 
tracker evolution by including the cosmic curvature K 
and test the viability of such a solution. In general the 
cosmological dynamics start from the regime away from 
the tracker, depending on the model parameters and ini- 
tial conditions. We shall also study such general cases 
and search for the model parameter space consistent with 
observational constraints. 



II. GALILEON COSMOLOGY 

The covariant Galileon gravity is described by the ac- 
tion m 



s = 



+ d 4 xC M , (1) 



where g is a determinant of the metric tensor g^ Vl M p \ 
is the reduced Planck mass, and Cj are constants. The 
covariant Lagrangians Ci (i = 1, ■ ■ ■ ,5) that respect the 
Galilean symmetry in the limit of the Minkowski space- 



time are given by 

Ci = M 3 <f>, C 2 = (V0) 2 , C 3 = (□^)(V^)7M 3 , 
£4 = (V</>) 2 [2(D</>) 2 - 2<p.^ v - i?(V</>) 2 /2] jU\ 
C 5 = (V</>) 2 [(D</>) 3 - 3(D0) <j>.^ v 

+2^/^^;/ - H,^ v ^ p G vp \lM <A , (2) 

where a semicolon represents a covariant derivative, M 
is a constant having a dimension of mass, and G vp is the 
Einstein tensor. For the matter Lagrangian Cm we take 
into account perfect fluids of radiation (energy density p r , 
equation of state w r = 1/3) and non-relativistic matter 
(energy density p m , equation of state w m = 0). 

We consider the FLRW space-time with the line ele- 
ment 



ds 2 = -dt 2 + a 2 (t) 



dr 2 



1-Kr 2 



\d9 2 



sm 



(3) 

where a(t) is the scale factor with the cosmic time t. The 
closed, flat, and open geometries correspond to K > 0, 
K = 0, and K < 0, respectively. Variation of the action 
([1]) with respect to g pv leads to the following equations 
of motion 



ZM'^H 2 



PDE + Pm + Pr + PK , (4) 

3M 2 ,H 2 + 2M 2 l H = -P DE - p r /S + p K /3 , (5) 



h3i?p„ 
AH Pr 



0. 



0. 



(6) 
(7) 



where px = SKAI^/a 2 , a dot represents a derivative 
with respect to t, H = a/a is the Hubble parameter, and 

Pde = - Cl M 3 (j)/2 - c 2 2 /2 + 3c 3 H<j> 3 /M 3 

-45c4H 2 <j> A /(2M 6 ) + 2lc 5 H 3 <j> 5 /M 9 , (8) 
'de = ciM 3 <j)/2 - c 2 tf/2 - c 3 ^ 2 4>/M 3 



+3c i( j) 3 [8H^ + (3H 2 + 2H)<j)]/(2M 6 ) 
-3c 5 i7^ 4 [5iJ0 + 2(H 2 + H)<j)]/M 9 . (9) 

Note that Eqs. © and are the generalization of 
those derived in Ref. [23| with an account of the cosmic 
curvature K. The closed-form equations for (f> and H 
can be derived by taking a time-derivative of Eq. Q and 
by combining it with Eq. ([5]). Since we are interested 
in the case where the late-time cosmic acceleration is re- 
alized by the field kinetic energy, we set c\ = in the 
following discussion. In this case the only solution in the 
Minkowski background {H = 0) corresponds to cj) = for 

The dS solution (H = H^s = constant) can be present 
for cb = 0ds = constant. We normalize the mass M to 
be M 3 = M p iH 2 s , which gives M w Kr 40 Af pl for H ds w 



10- bu M pl . Defining x dS 



"as 



/(ffdsMpi), Eqs. O and 



give the following relations at the dS point: 

c 2 x 2 ds = 6 + 9a - 12/3 , 
c 3 xl s = 2 + 9a - 9/3 , 



(10) 
(11) 
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where 



_ 4 



13 = c 5 x 5 dS . 



(12) 



The use of a and /3 is convenient because the coefficients 
of physical quantities and dynamical equations can be 
expressed by those variables. We note that the relations 
(fT0|) and (fTTj) are not subject to change under the rescal- 
ing Xds — ► 7^ds an d c; — > Ci/^f, where 7 is a real con- 
stant. Hence the rescaled choice of c,; can provide the 
same physics. If we use the parameters a and /?, such 
apparent different cases can be treated in a unified way. 
In Refs. [23|,[24[ the authors derived the viable parameter 
space on the (a, (3) plane in which the conditions for the 
avoidance of ghosts and Laplacian instabilities of scalar 
and tensor perturbations are satisfied for K = 0. 

In order to study the cosmological dynamics, it is use- 
ful to define the following dimensionless variables: 



'•1 



6dS#dS 

4>H 



>2 




(13) 



At the dS solution one has n = 1 and r 2 = 1. We define 
the dark energy density parameter 



DE 



Pde 



SM^H 2 



-(2 + 3a - 4/3)r?r 2 /2 + (2 + 9a - 9/3)r 2 r 2 
-15anr 2 /2 + 7/3r 2 , (14) 



where we have used Eqs. (fT0|l and (fTTj) . Then Eq. (j4j) can 
be written as 



f^DE + + + "jr = 1 



(15) 



where Q m = p m /(3M 2 l H 2 ), fi r = p r / '(3M^H 2 ) , and 
n K = P k/(3M^H 2 ) = —Kj (aH) 2 . 

From Eqs. (HI)-® we obtain the following autonomous 
equations for the variables r\, r 2 , fi r , and fix- 



K 



A 



(n - 1) n [n (n(-3a + 4/3 - 2) + 6a - 5/3) - 5/3] 



x [2 (fi r - ft K + 9) + 3r 2 (r?(-3a + 4/3 - 2) + 2r 2 (9a - 9/3 + 2) - 15^0 + 14/3)] , 



(16) 



1 



-[r 2 (6rf (r 2 (45a 2 - 4(9a + 2)/3 + 36/3 2 ) - (O r - O x - 7)(9a - 9/3 + 2)) + r?(-2(fi r - ft K + 33) 



x (3a - 4/3 + 2) - 3r 2 (-2(201a + 89)/3 + 15a(9a + 2) + 356/3 2 )) - 3na(-28fi r + 28fl K + 123r 2 /3 + 36) 
+10j8(-lin r + llfl K + 21r 2 /3 - 3) + 3r^r 2 (9a 2 - 30a(4/3 + 1) + 2(2 - 9/3) 2 ) + 3rfr 2 (3a - 4/3 + 2) 2 
+3r?r 2 (9a-9/3 + 2)(3a;-4/3 + 2))], (17) 

~H 

IT 



-4-2- 
-2-2- 



(18) 
(19) 



where a prime represents a derivative with respect to N = In a, and 



A = 2rfr 2 [72a 2 + 30a(l - 5/3) + (2 - 9/3) 2 ] + 4rf[9r 2 (5a 2 + 9a/3 + (2 - 9/3)£) + 2(9a - 9/3 + 2)] 

+4rJ[-3r 2 (-2(15a + l)/3 + 3a(9a + 2) + 4/3 2 ) - 3a + 4,3 - 2] - 24r!a(16r 2 /3 + 3) + 10/3(21r 2 /3 + 8) . (20) 
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The Hubble parameter follows from the equation 



IT 



5_r[ 

in 



4r 2 



(21) 



The solutions to Eqs. (|T5|) and (fT9"| are given 

by n r (N) = n [ r ] e- AN H 2 /H 2 (N) and Q K (N) = 

e~ 2N H 2 / H 2 (N) respectively, where the subscript 
"(0)" represents the values today (N — and a = 1). 
Hence these variables are related with each other via 



n K (N) = nf(n r (N)/n^)e 2N . 



(22) 



A. Tracker solution 



Equation (fT6"|) shows that there is an equilibrium point 
characterized by 



along which 



n = l 



^DE = ^2 , 



(23) 



(24) 



Originally the existence of the tracker solution (|23|) was 
found in Refs. [HHI for the fiat Universe (K = 0). The 
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above results show that the tracker is also present for 
K ^ 0. 

The epoch at which the solutions reach the tracking 
regime r± ~ 1 depends on model parameters and initial 
conditions. The approach to this regime occurs later for 
smaller initial values of r\ . If r\ < 2 initially, numerical 
simulations show that the solutions approach the tracker 
with the late-time cosmic acceleration. Meanwhile, for 
the initial conditions with r\ > 2, the dominant contri- 
bution to 57de comes from the Lagrangian £ 2 , so that 
the field energy density decreases rapidly as in the stan- 
dard massless scalar field. In the latter case the solutions 
do not get out of the matter era that starts from the 
radiation-matter equality [13, HH ■ 

From Eqs. (fT7 l) - (fT9| the variables r 2 , fl r , and fl x sat- 
isfy the following equations along the tracker: 



ft 



"if 



2r 2 (3 - 


3r 2 + fl r - 


fl K ) 


(25) 




l + r 2 






- fl K ~ 1 - 


7r 2 ) 


(26) 




1 + r 2 


* 


fix {flr 


- fl K + 1 - 


-5r a ) 


(27) 




l + r 2 





These equations do not have any dependence on a and 
(3. Combining Eqs. (f25|) and (f26|) . we obtain 
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8- 



2fl' r 

fi r 



Integrating this equation, it follows that 

r 2 = d ia 8 fl 2 r , 

where d\ is a constant. 

From Eqs. and (J2ZJ) we have 

fix fir 



(28) 



(29) 



(30) 



which is integrated to give 



f}K_ 

fir 



d 2 a 2 



with 



fl (0) 

d, = 



fl) 



(0) 



(31) 



Substituting Eqs. (gSJ) and into Eq. ([H, w e ob- 
tain the following integrated solution 



-1 + da a 



-l + d 3 a-<i 2 a 2 ) 2 



2dia 8 



(32) 

where d 3 is another constant. Note that another solution 
of fl r (i.e. minus sign in front of the square root) is 
not cosmologically viable, because of the divergence of 
fl r as a — > 0. The density parameter (|32p evolves as 
fl r ~ 1 + c^a in the early time (a <C 1). This demands 
the condition d 3 < provided that flr>E > 0. 



Using the density parameters fl m , fl 



(0) o (0) 



and fix to- 



day, the constants d\ and d% can be expressed as 



di = 



1 - Q 



(0) 



fl 



(0) 
K 



<fo = — 



(0) 



(33) 



(flr 0) ) 2 fl r 0) 

where we have used Eqs. ([T5]l. (HH , and (|3"2"]1. In the high- 
redshift regime (2> 1) the density parameters behave as 
fl DE ~ di/[(l+z) 6 {l+z-da)% fl r ~ (l + z)/(l + z-d 3 ), 
and fi^ ~ d 2 /[(l + z)(l + z - d 3 )]. 

Since i7 r oc p r /H 2 oc l/(a 4 H 2 ), it follows that 
H 2 /Hi = (fl r 0) /fl r )(l/a 4 ). Using Eqs. fl32) and (J33J), 
the Hubble parameter can be expressed in terms of the 
function of the redshift z = 1 /a — 1 : 



H(z) 

Hn 



in£>(i + z) 2 + 1qW(i + z) 3 + in(«»(i + z) 4 



i-n£>-n< > 



(0) 



fig ) +^(l + z) + ^(l + z) 



,(0) 



(34) 



This analytic estimation is useful to constrain the tracker 
solution from a number of observations. 



On the tracker the equation of state of dark energy 
wde = Pde/pde and the effective equation of state 



w c s = — 1 — 2H/(3H 2 ) arc given by 



W c ff 



fl r — fix + 6 
3(r 2 + 1) 
fl r — fix ~ 6r 2 
3(r 2 + 1) 



(35) 
(36) 



During the cosmological sequence of radiation (fl r ~ 1, 
\flK\ < 1, r 2 < 1), matter (O r < 1, \fl K \ < 1, 
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log 10 (l+z) log 10 (l+z) 

Figure 1: The evolution of the dark energy equation of state wde versus the redshift z for 0,^ = —0.1 (left) and S ]^ = 0.1 
(right) with a = 0.3 and /3 = 0.14. We choose the initial conditions as those given in the caption of figure 3 of Ref. [24j. The 
tracker solution (ri = 1) is shown as a bold line. For given Q, r (N), the curvature density parameter £Ik{N) is determined 
according to Eq. (|22[) . The two vertical lines denote the present time (z = 0) and the epoch at decoupling z = z», while the 
horizontal line corresponds to the ACDM model (u>de = — !)• 




-10 12 3 4 

log 10 (l+z) 



Figure 2: The evolution of woe for a — 0.3 and /3 = 
0.14 with various values of the cosmic curvature: = 
(—0.25,-0.1,0,0.1,0.25). The initial conditions are chosen 
to be ri = 1.5 x 10" 10 ,r 2 = 2.667 x 10~ 12 ,a- = 0.999992 at 
z = 3.63 x 10 8 [where N = - ln(l + z)]. The meaning of two 
vertical lines and the horizontal line is the same as in Fig. [1] 

r 2 < 1), and dS (O r < 1, \Cl K \ < 1, r 2 = 1) eras, 
the dark energy equation of state shows peculiar evolu- 
tion: iutje = —7/3 — > —2 — > —1, whereas the effective 



equation of state evolves as w e g = 1/3— » — > — 1. If 
the curvature density parameter Qk is non-negligible to- 
day, this also gives some contribution to u>de and w; e ff 
during the evolution from the end of the matter era to 
the present epoch. In Fig. [T] the evolution of the tracker 
solution (shown as a bold line) for 0^ = —0.1 does not 

look very much different from that for Qr^' =0.1. We 
caution, however, that inclusion of the cosmic curvature 
changes the diameter distance as well as the luminosity 
distance relative to the flat Universe. 

The conditions for the avoidance of gho sts and insta- 
bilities have been derived in Refs. [23l. |24| for D,k = 0. 
Even in the presence of the cosmic curvature it is a good 
approximation to use the results in the flat case, because 
fix is much smaller than 1 in most of the expansion his- 
tory of the Universe. Hence we shall use the allowed re- 
gion in the (a, /3) plane shown in the figure 1 of Ref. [23| . 



B. General solutions 

There is another case in which the solutions start to 
evolve from the regime ri <C 1. If ri <C 1 initially, the 
ghosts are absent for j3 > [H,^. Provided that n <C 1 
the propagation speeds of scalar and tensor perturbations 
are positive during radiation and matter eras, so that no 
instabilities are present. 

In the regime n < 1, the variables n and r 2 satisfy 
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the following approximate equations 

, 9 + Q, r - n K + 21/3r 2 

8 + 21/?r 2 ri ' 
, 3 + llfi r - l\n K - 21/3r 2 
fa ~ 8 + 21/3r 2 



-T2 ■ 



(37) 
(38) 



As long as {/3r 2 , <iC 1, the evolution of r\ and r 2 

during the radiation (matter) era is given by r\ oc a 5 / 4 
and r 2 oc a 7 / 4 (ri ex a 9 / 8 and r 2 c>c a 3 / 8 ). Hence the field 
velocity evolves as (f> cx i 3 / 8 during the radiation era and 
4> oc t 1 / 4 during the matter era. Note that the evolution 
of <j> is slower than that for the tracker (i.e. </> oc t). 

In Rcf. [24[ it was shown that the tracker is stable in 
the direction of t\ by considering a homogeneous pertur- 
bation Sri . This means that once the solutions reach the 
tracker the variable r\ does not repel away from 1. The 
epoch at which the solutions approach the tracker regime 
depends on the initial values of T\. 

The dark energy equation of state u>de and the ef- 
fective equation of state w e s m the regime n *C 1 arc 
approximately given by 



1"DE 



1 + Q r - £1 



K 



8£l r 



±21/3r 2 
8n K - 21/3r 2 



3(8 + 21/3r 2 ) 



(39) 
(40) 



Provided that {(3r 2 , \£Ik\} <C 1, one has wde — 

— l/4,w e g ~ 1/3 during the radiation era and wde — 

— l/8,w e ff — during the matter era. The evolution of 
u>de in the regime r\ <C 1 is quite different from that for 
the tracker solution. 

In Fig. Q] we plot the evolution of u>de for a number 
of different initial conditions. Since n < 1 initially, the 
solutions start to evolve from the value iude — —1/4 in 
the radiation era. For larger initial values of r± they ap- 
proach the tracker earlier. This tracking behavior occurs 
irrespectively of the signs of In Fig. [5] we find that 
the effect of the cosmic curvature slightly modifies the 
evolution of u>de in the low-redshift regime. 



III. METHOD OF LIKELIHOOD ANALYSIS 

In this section we show the method of our likeli- 
hood analysis to place observational constraints on the 
Galileon cosmology discussed above. The modified back- 
ground cosmological evolution in this theory affects the 
diameter distance to the last scattering surface as well 
as the luminosity distance. The modification from the 
ACDM model can be tested by using the data of CMB, 
BAO and SN la. 



A. CMB shift parameters 

The positions of CMB acoustic peaks are affected by 
the expansion history of the Universe from the decoupling 



epoch to today. In order to quantify the shift of acous- 
tic peaks we use the data points (l a ,lZ,z*) of Ref. || 
(WMAP_7_h where l a and K are two CMB shift parame- 
ters [4(H49| ] and z* is the redshift at decoupling. For the 
FLRW metric © we have [H 



K 



(0) 



sinh 



fi 



dz 



H(z)/H c 



(41) 



Numerically it is convenient to integrate the following 
equation from the redshift z = 0toz = z t : 

an _ i 

d7 ~~ H(z)/H Q 



n 



(0) 



n 2 n 



(0) 



K 



(42) 



The multipole l a is defined by l a = ■Kda{z*)/r s {z >i 

where Sa\z*) 
ing angular diameter 
J z °° dz/[^/3(l + R s (z))H(z)] is the sound horizon at the 



distance 



is 
and 



the comov- 

r s {z*) = 



decoupling, 
where tt^ and tt^ 



Note that R s (z) = 3ft[ 0) /[4f^ 0) (1 + z )], 

are the today's density parameters 
of baryons and photons respectively. We neglect the 
contribution of dark energy and the cosmic curvature 
for z > z* to estimate r s (z*). The dark energy density 
parameter can be in fact neglected even for non-tracker 
solutions with ri < 1, because Ode — 7/3r 2 decreases 
toward the past (r 2 oc a 3 / 8 during the matter era and 
r 2 oc a 7 / 4 during the radiation era). It then follows that 
(see e.g., [U) 

-l 



In 



3tt 



y/Rs{z,) + R s (z cq ) + y/T+ fl a (z.) 
l + ^i? s (z eq ) 



(0) 



K. 



(43) 



For the redshift z* there is a fitting formula by Hu and 
Sugiyama [5^ ]: 

-0.738\ 



1048 1 + 0.00124^: 



(44) 



where g\ 
0.560/ (1 H 



0.0783w b - a238 / (1 



39.5u;0- 763 ), 52 



21.1wJ- 81 1 



LJb = Q^h 2 , and uj„ 
(h correspond to the uncertainty of the Hubble param- 
eter H today, i.e. H = 100 h km sec" 1 Mpc" 1 ). The 
redshift at the radiation-matter equality is given by 



/ 



1. 



For a flat prior, the 7-year WMAP data (WMAP7) 
measured best-fit values are [f| 



CMB — 









■I 







302.09 ±0.76 



1091.3 ±0.91 
The corresponding inverse covariance matrix is Q 



r 1 



2.305 29.698 -1.333 
29.698 6825.270 -113.180 
-1.333 -113.180 3.414 



(45) 



(46) 
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Wc thus define 



L CMB 



l a - 302.09 
U - 1.725 
z* - 1091.3 



where cIa{z) is the diameter distance defined by 
1 1 . , f fZim f z dz 



(47) 



d A {z) 



1 + z H Jnf 



: sinh 



o H(z)/H 



and construct the contribution of CMB to x 2 as 

XCMB = ^CMB CcMB ^CMB • (48) 

Notice that Xcmb depends on the parameters (flm , 
fi£ , h). In the analysis of the Galileon model we use the 
priors h = 0.71 and n { ° ] = 0.02258 h~ 2 @. As we shall 
see later on, varying h does not affect the final results of 
our study. The density parameter of radiation today is 



(53) 

This encodes the visual distortion of a spherical ob- 
ject due to the non-euclidianity of the FLRW space-time, 
which is equivalent to the geometric mean of the distor- 
tion along the line of sight and two orthogonal directions. 

We thus construct 



XbAO 



^M)~ 0.1905' 
TJTfefe- 0.1097 ^ 



and using the inverse covariancc matrix 



= fi( >(l + 0.2271JVeff), 



(49) 



u BAO 



where is the photon density parameter and N e g is 

the relativistic degrees of freedom. We adopt the stan- The contribution of BAO to \ 2 is 



30124 -17227 
-17227 86977 



dard values Q 



(0) 



5]. 



2.469 x 10" 



and N eB = 3.04 



Xbao 



— -^BAO^RAO-^BAO 



(54) 



(55) 



(56) 



In Ref. [46j Wang and Mukherjee placed observational 
constraints on 1Z and l a for several different dark en- 
ergy models: the ACDM model, constant wde mod- 
els, and the models described by the parametrization 
ujfje = wo + w a (l — a). They showed that the resulting 
bounds on 1Z and l a are independent of the dark energy 
models. While our Galileon model does not belong to 
the models described above, the method using the two 
distance measures 1Z and l a is expected to be trustable 
as well. In fact the Galileon model is consistent with 
a number of assumptions [53j that validate the analysis 
using the two WMAP distance priors. 



B. BAO 

For the BAO wc apply the maximum likelihood 
method [H] using the data points of Ref. Q (SDSS7): 



C. SN la 

The analysis of SN la standard candles is based on the 
method described in Ref. (48|. We will mainly use the 
Constitution SN la dataset of Hicken et al. [55|], which 
constitutes in total of 397 SN la. We will also use the 
recently released update to the Union set, i.e. the Union2 
dataset (56|. 

The SN la observations use light curve fitters to pro- 
vide the apparent magnitude m(z) of the supernovae at 
peak brightness. This is related with the luminosity 
distance di,(z) through m(z) = M + 5 log 10 (di,/10 pc), 
where M is the absolute magnitude. Note that the lumi- 
nosity distance is given by 



d L (z) = {l + z) 2 d A {z), 



(57) 



V 



BAO 



Dy (0.2) 

r 3 (z d ) 
Dv(0.35) 



0.1905 ±0.0061 
0.1097 ±0.0036 



(50) 



where r s (zd) is the sound horizon at the baryon drag 
epoch Ztf. For z^ we use the fitting formula by Eisenstein 
and Hu (H: 



1291w 



1 + 0.659o;0: 828 



(l + byuf) , (51) 



where dA{z) is the angular diameter distance defined in 
Eq. (|53p . Defining the dimcnsionlcss luminosity distance 
as <1l(z) = / Hq 1 , the theoretical value of the ap- 

parent magnitude is 

m th (z) = M(M, H ) + 5log 10 {d L (zj) , (58) 

where M = M - 5 log 10 h + 42.38 [H]. 

The theoretical model parameters are determined by 
minimizing the quantity 



where h = 0.313 w m 419 (l + 0.607 w^ 674 ) and b 2 = 
0.238 uj° m 223 . The dilation scale D v at the redshift z is 



N 

xlNIa(^n° ) ,P J )=E 



[Hobs(Zi) - ^th(Zi)f 



(59) 



fl I 



D v (z) 



(l + z) 2 d 2 A (z) 



H(z) 



1/3 



(52) 



where N is the number of the SN la dataset, pj is the set 
of parameters to be fitted, and i are the errors due to 
flux uncertainties, intrinsic dispersion of SN la absolute 
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magnitude and peculiar velocity dispersion. These er- 
rors are assumed to be Gaussian and uncorrelated. The 
theoretical distance modulus is defined as 

Mth(z») = rn th (z t ) - M = 5log 10 (d L (z)) + /j, , (60) 

where (j,q = 42.38 — 51og 10 /i. The steps we have fol- 
lowed for the usual minimization of (|59[) in terms of its 
parameters are described in detail in Refs. [5714591 ]. 

D. Two information criteria 

In order to see whether the Galileon model is favored 
over the ACDM model, we will also use the two infor- 
mation criteria known as AIC (Akaike Information Cri- 
terion) and BIC (Bayesian Information Criterion), see 
Ref. [60] and references there in. The AIC is defined as 

AIC = -2 In £ max + 2k , (61) 

where the likelihood is defined as C oc e~ x ' , the term 
—2 ln£ max corresponds to the minimum % 2 , and k is the 
number of parameters of the model. The BIC is defined 
similarly as 

BIC = -21n£ max + fclnAT, (62) 

where N is the number of data points in the set under 
consideration. 

According to these criteria a model with the smaller 
AIC/BIC is considered to be the best and specifically, 
for the BIC a difference of 2 is considered as positive 
evidence, while 6 or more as strong evidence in favor of 
the model with the smaller value. Similarly, for the AIC 
a difference in the range between and 2 means that the 
two models have about the same support from the data 
as the best one, for a difference in the range between 2 
and 4 this support is considerably less for the model with 
the larger AIC, while for a difference > 10 the model with 
the larger AIC is practically irrelevant [6(J Hl[ . 

IV. OBSERVATIONAL CONSTRAINTS 

In this section we present the observational constraints 
on the Galileon cosmology. We first consider the tracker 
solution and later we study general solutions discussed 
in Sec. Ill Bl In the latter case we shall explore the whole 
parameter space in terms of (a, /?) constrained by the 
conditions for the avoidance of ghosts and instabilities 
(figure 1 in Ref. (H). 

A. The tracker 

As we showed in Eq. (|34|). the Hubble parameter for 
the tracker is known as a function of the redshift z. Since 
the value of fi[ 0) is fixed from the CMB [see Eq. (03)] , the 



0.02 




-0.02 

-0.03 ' — ■ ■ ■ 1 ■ ■ ■ 1 ■ ■ ■ 1 ■ ■ ■ 1 ■ — 

0.24 0.26 0.28 0.30 0.32 

Figure 3: The 68.3% (la) -95.4% (2a) \ 2 confidence contours 
in the (flm , Q,^ ) plane for (i) the tracker model and (ii) the 
ACDM model. Both contours correspond to the combination 
of all three datasets, i.e. Constitution SN Ia+BAO+CMB. 
The best fit parameters are found in Table [I] The difference 
in x 2 between the two models is S\ 2 ~ 22, which corresponds 
to ~ 4.3a. Hence the tracker model is severely disfavored 
with respect to the ACDM model. 

tracker has only two free parameters ft^ and 0,^ ■ I n 
what follows the priors h = 0.71 and fij; 0) = 0.02258 /i" 2 
are used in order to simplify the analysis. We have 
checked that having ft, as a free parameter for the likeli- 
hood analysis of the CMB and BAO does not change the 
results much, as the best fits are always within 1 %. 

We should also mention that we have not solved the 
full perturbation equations for the Galileon field in or- 
der to calculate the effect on the CMB. However, the 
modified growth of perturbations only affects the large- 
scale CMB spectrum (the multipoles / < 10) through the 
Integratcd-Sachs- Wolfe (ISW) effect. Usually this does 
not provide tighter constraints than the CMB distance 
measures explained in Sec. IIII Al Hence we expect that 
the results are not subject to change much. 

In TableUwe show a comparison of the tracker given in 
Eq. (|34|) to the ACDM model for various data combina- 
tions. If we only use either of the SN la data (Constitu- 
tion or Union2), the x 2 f° r the tracker is similar to that in 
the ACDM for both the flat and non-flat cases. Hence the 
two models are practically indistinguishable from each 
other. The tracker solution is also consistent with the 
individual observational constraint from either CMB or 
BAO. The combined data analysis of CMB+BAO shows 
that, as long as we take into account the cosmic curva- 
ture, the x 2 for the tracker is even smaller than that in 
the ACDM. 

In Fig. 02 we plot the X 2 (= x|n + Xcmb + Xbao) 
confidence contours in the (flm , f2j£ ) plane at the 
68.3% (la) - 95.4% (2a) levels. Both contours corre- 
spond to the combination of all three datasets: Consti- 
tution SN Ia+BAO+CMB. As we see in Table U the 
difference in x 2 between the two models is 5x 2 ~ 22. 
This corresponds to ~ 4.3er, and the tracker solution is 



9 



Table I: Comparison of the tracker solution characterized by (|34|) to the ACDM model for various data combinations. If we use 
the SN la data alone, the model is practically indistinguishable from the ACDM. However, if we add the CMB+BAO data, 
the tracker solution is disfavored relative to the ACDM. In all cases the error estimation is only statistical. In the cases where 
— it means that we assume the flat Universe in the data analysis. 



Model (datasets) 


Xmin 


n (0) 


fi (0) 


ACDM (CMB) 


0.22 





i\ r\ f en i r\ aao 

0.268 ± 0.003 


Tracker (CMB) 


16.50 


() 


0.316 ± 0.004 


ACDM (CMB) 


0.21 


0.000 ± 0.004 


0.27 + 0.01 


Tracker (CMB) 


0.21 


-0.016 + 0.003 


0.27 + 0.01 


ACDM (BAO) 


1.58 





r\ onp + 0.047 

0.306I 041 


Tracker fBA(Yl 


0.47 


o 


263+ 039 

U.ZUO_q 034 


ACDM (BAO) 


0.00 


-0.37 + 0.23 


0.22 + 0.06 


Tracker (BAO) 


0.00 


-0.33 + 0.61 


0.22 + 0.07 


ACDM (CMB+BAO) 


2.65 





0.268 ± 0.004 


Tracker fCMB+BAOl 


18.71 


o 


316 + 004 


ACDM (CMB+BAO) 


2.56 


0.001 + 0.004 


0.27 + 0.01 


Tracker (CMB+BAO) 


0.67 


-0.016 + 0.26 


0.27 + 0.12 


ACDM (Constitution) 


465.512 





0.289 ± 0.022 


TVapVpr i f^nn tit i tnt inn l 


465 508 


o 


392 + 024 


ACDM (Constitution) 


465.353 


-0.09 + 0.23 


0.325 + 0.092 


Tracker (Constitution) 


465.423 


0.076 ± 0.258 


0.358 + 0.118 


ACDM (Umon2) 


541.011 





0.269 ± 0.020 


TV^ pkpt* ( T Tti i r>T~i 1 


541 088 


o 


369 + 099 


ACDM (Union2) 


540.879 


-0.07 + 0.19 


0.296 ± 0.076 


Tracker (Union2) 


540.853 


0.106 + 0.217 


0.324 + 0.097 


ACDM (Constitution+BAO) 


467.214 





0.269 ± 0.020 


Tracker ('nnnstitntinn^R AO^l 


472.729 


o 


367 + 099 


ACDM (Constitution+BAO) 


466.721 


-0.052 ± 0.074 


0.306 ± 0.028 


Tracker (Constitution+BAO) 


467.066 


0.172 + 0.067 


0.311 + 0.029 


ACDM (Umon2+BAO) 


543.210 





0.277 ± 0.018 


Tracker ('TIninn24-R AO 1 ) 


546 499 


o 


349 + 090 


ACDM (Union2+BAO) 


542.013 


-0.077 + 0.071 


0.296 + 0.027 


Tracker (Union2+BAO) 


542.182 


0.148 + 0.066 


0.303 ± 0.028 


ACDM (Constitution+BAO+CMB) 


469.024 





0.269 ± 0.004 


Tracker (Constitution+BAO+CMB) 


494.397 


o 


318 ± 004 


ACDM (Constitution+BAO+CMB) 


468.543 


0.002 + 0.003 


0.275 + 0.010 


Tracker (Constitution+BAO+CMB) 


490.225 


-0.008 ± 0.003 


0.296 + 0.011 


ACDM (Union2+BAO+CMB) 


543.660 





0.268 + 0.003 


Tracker (Union2+BAO+CMB) 


565.576 





0.318 + 0.004 


ACDM (Union2+BAO+CMB) 


543.582 


0.001 + 0.003 


0.270 + 0.010 


Tracker (Union2+BAO+CMB) 


560.183 


-0.008 ± 0.003 


0.294 + 0.011 
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severely disfavored with respect to the ACDM. A simi- 
lar conclusion is reached from the combined data analy- 
sis of Union2-r-BAO-|-CMB, which gives the difference of 
S\ 2 ~ 16 relative to the ACDM model. 

The reason why the tracker is disfavored can be ex- 
plained by inspecting Table U carefully. For the flat Uni- 
verse the CMB data are practically incompatible with 
the tracker because of the peculiar evolution of i«de- 
The main effect on the modification of the angular di- 
ameter distance to the last scattering surface comes from 
the contribution in the low-redshift regime (z < a few). 
In the non-flat background the CMB data are consistent 



with the tracker with the best fit values ft 



(0) 



0.27 and 



ftp ~ —0.016. This value of ftm is quite smaller than 



the SN la best fit > 0.32), with either the Consti- 

tution or the Union2 datasets. 

Regarding the SN la observations, the maximum red- 
shift for the Constitution dataset is z max = 1.551 and 
z max = 1.4 for the Union2 dataset, while the average 
rcdshifts are z mean = 0.4 and z mea ,n = 0.35 respectively. 
At z = 0.4 and z = 1.4 the dark energy equations of state 
are given by ujde ~ —1.4 and wde ~ —1.9, respectively. 
Since these values are away from u>de = — 1, the tracker 
can be compatible with the SN la data at the expense of 
choosing the values of ftS larger than 0.3. 

For the BAO datasets the best-fit value of Si™ in the 
non-flat Universe is even smaller (ftm ~ 0.22) than that 
constrained by the CMB shift parameters. Thus, as each 
dataset favors the value of fl^m in a different range, the 
combined analysis with all datasets does not favor the 
tracker solution. 



B. General solutions 

Next, we proceed to observational constraints on the 
general solutions to Eqs. (|T6]) - ([T9^) . Unlike the tracker 
solution, this case depends on the parameters a and /3 
as well as the initial conditions of ri, ra, and fl r . Note 
that the initial condition of flK is fixed by using Eq. (1221) . 
From Eq. (|2~lj) the Hubble parameter H(N) is integrated 
to give 



H(N) /n(0) 



ri(N) 



5/4 



7-2(0) 

r 2 (N) 



1/4 



(63) 



which can be used to compare the model with the data. 

We place observational bound on (a, {$) constrained 
by the conditions for the avoidance of ghosts and insta- 
bilities. We also restrict j3 > to avoid the appearance 
of ghosts at the early cosmological epoch [23|, |2J| . Since 
the tracker is disfavored by the data, the cosmological 
evolution in which the solutions approach the tracker at 
late times (i.e. smaller initial values of 7*1) is in general 
favored. For given a and /3, we search for the viable ini- 
tial conditions of n, r 2 , and fi r that can be consistent 
with the data. We then carry out the likelihood analysis 
by varying a and ft to find the viable parameter space. 



p 




Figure 4: The evolution of the dark energy equation of state 
wde for the tracker solution (|34[) and the non-tracker case of 
Table [TT] [denoted by an (*), i.e. a = 1.862 and = 0.607 
with the initial conditions given in the text]. The meaning 
of two vertical lines and the horizontal line is the same as in 
Fig. [I] 



In Table UU wc show the best fit and the best fit pa- 
rameters along their la errors and the AIC/BIC for sev- 
eral fixed values of ftjp (= 0, —0.01,0.01). According to 
the AIC statistics, the ACDM and the non-tracker solu- 

0.01 



(0) 

K 



tion with the positive curvature parameter ft 
(with the best-fit values a = 1.862 and (3 = 0.607) have 
more or less the same support. For this set of best- 
fit parameters we can always find a set of initial con- 
ditions for r% and r 2 that lead to the late-time tracker 
allowed by observations. For example, for the param- 
eters a = 1.862 and (3 = 0.607, the initial conditions 
giving rise to the late-time tracker allowed by obser- 
vations are n = 1.50 x 10" 10 ± 3.88 x 10" 12 , r 2 = 
2.67 x 10" 12 ±2.17x 10~ 13 (95 % CL), for ft r = 0.999992 
at the redshift z = 3.63 x 10 8 . These conditions can 
be translated into the initial condition for <f> by using 
Eq. KB!}. In the case of a = 1.862 and /3 = 0.607 it follows 
that 4.38 x 10~ 6 < 0/0 d s < 4.57 x 10~ 6 at z = 3.63 x 10 8 , 
thus requiring a certain amount of fine tuning. 

For ftP = -0.01 the data support the ACDM model 
more, as the fit is particularly bad. On the other hand, 
according to the BIC statistics, there seems to be strong 
evidence in favor of the flat ACDM over all non-tracker 
cases as ABIC > 6. 

In Fig.[l]we illustrate the evolution of wde for the solu- 
tion denoted by an asterisk in Table [Til together with the 
tracker solution. The former corresponds to the model 
parameters a = 1.862 and j3 = 0.607 with the initial 



conditions n = 1.5 x 10 



-10 



'2 



2.667 x 10" 12 , and 
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Table II: Comparison of general solutions of Eqs. ()16|) - (|19[ ) to the ACDM model for various data combinations and values of 
the curvature parameter Qjp . For the model that is denoted by an asterisk (*), we show the evolution of iude in Fig. [4] (labeled 
as "non-tracker"). 



Model/Data 



Xmin 



Best fit parameters AIC AAIC AIC ABIC 



(Constitution+BAO+CMB) 





ACDM 


= 0) 


469.024 


fi^ = 0.269 ± 0.004 


471.024 





475.020 







ACDM (fi™ 




468.543 


n^ ) =0.002±0.003 
n^ , =0.267±0.004 


472.543 


1.519 


480.536 


5.515 




(a, /J), (ng> 


= 0) 


470.912 


a=1.411±0.056 
/3=0.422±0.022 


474.912 


3.888 


482.905 


7.884 




-0.01) 


537.543 


a=0.956±0.079 
,3=0.477±0.027 


541.543 


70.519 


549.536 


74.515 


(*: 


) (a, 13), (nP 


= 0.01) 


468.545 


a = 1.862±0.058 
/3=0.607±0.024 


472.545 


1.521 


480.538 


5.517 








468.311 


a=1.401±0.159 
,3=0.425±0.064 

f!^ ) =0.287±0.014 
a^'=-0.003±0.005 


476.311 


5.287 


492.297 


17.276 



(Union2+BAO+CMB) 



ACDM (fi^ = 0) 
ACDM (fi^ 5 / 0) 

(a,/3), =0) 
(a,/3), (fi^ = -0.01) 
(a,|8), (0^=0.01) 



d^ 7 " = 3 545.660 549.992 

547.582 1.922 556.245 6.254 

547.895 2.235 556.558 6.567 

614.919 69.259 623.582 73.591 

547.646 1.986 556.309 6.318 

550.290 4.630 567.616 17.625 



543.660 

543.582 

543.895 
610.919 
543.646 
542.290 



n^=o.ooi±o.oo3 

n^ ) =0.270±0.010 

a=1.404±0.057 
/3=0.419±0.023 

a=0.956±0.089 
,3=0.477±0.027 

a=1.8SS±0.058 
,3=0.603±0.024 



a=1.401±0.779 
,3=0.423±0.313 



n'° ) =0.282±0.016 



Co). 



-0.002±0.019 



fi r = 0.999992 at z = 3.63 x 10 s . In this case the so- 
lution does not reach the tracker by today and hence it 
does not exhibit strong phantom behavior. We find that 
the early tracking behavior is in fact disfavored by the 
combined data analysis. 

The results in Table |TT] for fixed fi)^ imply that the 
generic solutions fare quite well for the open Universe (i.e. 
positive fi^), whereas they are not in good agreement 

with the data for the closed Universe with 0,^ < —0.01. 
However, if we perform a general fit where all 4 pa- 
rameters (a, /?, fim , fij^ ) are allowed to vary, then a 
slightly negative value of fi^ is favored, e.g., fi^ = 
—0.003 ±0.005 for the Constitution set. The reason why 
the general solutions are not particularly favored in the 
AIC and BIC tests lies in the fact that the number of 
parameters is larger than those in the flat ACDM. There- 
fore, we conclude that the general solutions exhibit rich 
phenomenology that can in general be in good agreement 



with the combined observational constraints, while parts 
of the parameter space being in mild tension with the 
observations. 



V. CONCLUSIONS 

In this paper we have placed observational constraints 
on the covariant Galilcon cosmology by using the most re- 
cent data of SN la (Constitution and Union2 sets), CMB 
(WMAP7), and BAO (SDSS7). In this theory there is an 
interesting tracker solution that finally approaches a dc 
Sitter solution responsible for dark energy. By including 
the cosmic curvature fi]£ , we derived the analytic for- 
mula (p4|) about the evolution of the Hubble parameter 
for the tracker. This formula is convenient because we 
only need to vary the two density parameters fi^ an d 
fi^ for the likelihood analysis as in the ACDM model 
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(with the radiation density parameter £lf- fixed by the 
CMB). 

If we use either of the SN la data (Constitution or 
Union2) alone, the x 2 for the tracker is similar to that 
in the ACDM model. We also found that, as long as 
the cosmic curvature is taken into account, the tracker 
solution is compatible with the individual observational 
bound constrained from either CMB or BAO. However, 
the combined data analysis of Constitution+BAO+CMB 
shows that the difference of \ 2 between the tracker and 
the ACDM is Sx 2 ~ 22 (or ~ 4.3a). Hence the tracker 
is severely disfavored with respect to the ACDM. From 
the combined data analysis of Union2+BAO+CMB, we 
reached a similar conclusion: the difference of S\ 2 ~ 16 
relative to the ACDM. The reason for this incompatibility 
is that the SN la data favor the large values of Clm (> 
0.32), whereas the CMB and BAO data constrain smaller 

values of fi^ } (< 0.27). 

We also studied the general solutions to Eqs. (TTo]) - (rnj)) 
for the parameters a and /3 constrained by the condi- 
tions for the avoidance of ghosts and instabilities. In 
general the solutions that approach the tracker only at 
late times are favored from the combined data analysis. 
By choosing several fixed values of Cl^ , we found that 
the generic solutions can be consistent with the data for 
the open Universe (0^ > 0). For example, the general 

solutions with 0^ = 0.01 and the model parameters 
(a,/3) = (1.862,0.607) give the similar value of \ 2 to 
that in the the ACDM. In this case the AIC statistics 
also have the same support for the two models. For the 
models with largely negative £1% such as Q$ < —0.01 
the data favor considerably more the ACDM model, as 
the fit is particularly bad. 

The BIC statistics show that the general solutions, 
with all 4 parameters (a, /3, fim , ^k) are varied, are 
not particularly favored over the ACDM model (because 
ABIC > 6). This mainly comes from the statistical prop- 
erty that the number of model parameters is larger than 
those in the flat ACDM. In fact the general solutions 



with a non-zero curvature can be well consistent with 
the combined data analysis. 

It will be of interest to study the evolution of matter 
density perturbations to confront the Galileon cosmol- 
ogy with the observations of large scale structure (LSS). 
In particular, the presence of the field derivative cou- 
plings with R and G up in the expression of £4 and £5 
will change the effective gravitational coupling [62j . This 
should modify the growth rate of matter perturbations as 
well as the ISW effect on large-scale CMB anisotropics. 
Compared to the distance measures 1Z and l a we dis- 
cussed in this paper, the ISW effect is in general not so 
powerful to constrain dark energy models. However, in 
the case where the modification of the effective gravita- 
tional coupling is significant, the ISW effect may pro- 
vide a strong constraint. In addition, the LSS-ISW anti- 
correlation found in a similar Galileon-like model [28| can 
be a useful tool to constrain the parameter space (a, /3) 
further. We will leave this issue for a future work. 
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